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Abstract 


This  report  documents  the  third  stage  of  the  continuation  of  the  fully  coupled  numerical 
modeling  of  the  detonation  of  a  simplified  munitions  stack  in  a  temporary  storage  area  and  the 
subsequent  effects  on  the  immediate  surroundings  of  the  stack.  Three  plausible  configurations 
of  this  munitions  stack,  referred  to  as  the  “donor”  stack,  an  intervening  water  barricade,  and  an 
“acceptor”  munitions  stack,  are  modeled  in  two-dimensional  (2-D)  Cartesian  hydrocode 
computations  using  the  CTH  hydrodynamics  computer  code.  The  distance  between  each 
munitions  stack  and  the  barricade,  referred  to  here  as  the  “standoff”  distance,  is  varied  from  one 
computation  to  the  next,  with  the  physical  characteristics  of  the  munitions  stacks  and  barricade 
themselves  remaining  unchanged.  The  donor  stack  is  modeled  as  an  uncased,  condensed 
high-explosive  charge  with  a  rectangular  cross  section.  The  water  barricade  has  a  relatively  thick 
rectangular  cross  section,  and  the  acceptor  stack  is  modeled  as  a  solid  iron  rectangle.  The 
loadings  on  both  the  barricade  and  the  acceptor  stack  are  computed,  as  are  their  fully  coupled 
responses  to  those  loadings.  Only  a  relatively  weak  inverse  functional  relationship  with  standoff 
distance  was  found  in  the  barricade  response.  Weak  correlations  with  both  standoff  distance  and 
face  separation  were  also  found  for  all  parameters  that  were  evaluated  for  the  acceptor  stack 
response.  The  results  are  also  compared  with  those  of  the  first  two  parts  of  this  study  on  the 
coupled  blast  loading  and  response  computations  for  a  massive  water  barricade  with  a  trapezoidal 
cross  section  and  computations  for  a  thin  rectangular  water  barricade. 
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1.  INTRODUCTION 


This  report  documents  the  continuation  of  a  study  of  the  detonation  of  a  single  munitions 
stack  within  a  postulated  munitions  temporary  storage  area  and  the  subsequent  effects  on 
its  surroundings.  The  terminology  for  the  main  features  that  are  modeled  in  the  storage  area 
is  the  same  as  before.^’  ^  The  detonating  munitions  stack  is  referred  to  as  the  “donor”  stack; 
the  remaining  munitions  stacks  that  are  in  the  storage  area  and  subject  to  blast  loading 
from  the  donor  stack  are  the  “acceptor”  stacks;  a  postulated,  field-expedient,  protective  wall 
between  any  two  munitions  stacks  is  the  “barricade;”  and  the  distance  from  the  base  of 
a  munitions  stack  to  the  base  of  a  barricade  is  the  “standoff”  distance.  As  stated  before, 
the  primary  purpose  of  protective  barricades  is  to  prevent  a  direct,  line-of-sight  path  from 
existing  for  either  blast  or  fragments  between  munitions  stacks  in  proximity  to  one  another. 
Additionally,  the  impact  of  any  part  of  a  barricade  on  an  acceptor  stack  must  not  itself 
be  capable  of  initiating  an  exothermic  reaction  in  the  acceptor  stack.  This  computational 
study  is  one  part  of  a  larger  study  by  the  U.S.  Army  Research  Laboratory  (ARL)  on  behalf 
of  its  customer,  the  U.S.  Army  Defense  Ammunition  Logistics  (Ammolog)  Activity,  of  the 
dynamics  of  a  detonating  munitions  stack  and  the  effectiveness  of  field-expedient  barricades 
in  preventing  a  subsequent  chain  reaction  among  acceptor  stacks.  The  logic  for  considering 
the  development  of  field-expedient  barricades  has  been  previously  discussed  in  detail.  ^ 

The  first  computations  in  this  overall  study  were  two  independent,  “uncoupled”  compu¬ 
tations,  the  results  of  which  have  already  been  reported.^  The  first  of  those  of  computations 
modeled  only  the  detonation  of  the  donor  stack  and  the  subsequent  blast  loading  on  and 
response  of  a  massive  trapezoidal  water  barricade.  That  computation  was  run  until  the  barri¬ 
cade  had  achieved  a  nearly  steady  bulk  velocity  toward  the  location  of  an  acceptor.  However, 
that  computation  failed  to  continue  beyond  8  ms  simulated  time  because  of  numerical  sta¬ 
bility  problems.  Therefore,  it  was  necessary  to  have  a  second  computation  that  modeled 
a  reconstituted  barricade,  put  back  into  its  original  trapezoidal  cross-sectional  shape,  and 
traveling  toward  an  acceptor  stack  at  its  late-time  velocity  from  the  first  computation.  It 
was  started  at  the  instant  of  impact  of  the  barricade  on  the  acceptor  stack.  Thus,  these 
computations  are  termed  uncoupled  because  of  that  use  of  two  independent  computations. 
This  splitting  of  the  computational  problem  was  necessary  because  of  hydrocode  stability 
problems  in  the  very  difficult  breakup  phase  of  the  barricade  with  an  earlier  version  of  the 
hydrocode.  The  hydrocode  is  discussed  later  in  this  report.  That  first  uncoupled  study 
showed  that  the  massive  trapezoidal  water  barricade  was  effective  in  keeping  air  blast  load¬ 
ing  on  the  acceptor  stack  low  and  explosive  products  from  impinging  on  the  acceptor  stack. 
The  peak  pressures  on  the  acceptor  stack  face  caused  by  the  impact  on  the  reconstituted 
trapezoidal  water  barricade  were  approximately  500  MPa  (5  kbar).  As  with  any  uncoupled 
computations,  there  remained  open  questions  as  to  the  validity  of  performing  the  compu¬ 
tations  in  this  way  and  what  errors  might  have  been  introduced  because  of  that  approach. 
Also,  the  effects  of  variations  in  standoff  distance  were  yet  to  be  addressed. 
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After  this  initial  pair  of  computations,  three  separate  series  of  “fully  coupled”  compu¬ 
tations  for  water  barricades  were  performed.  These  computations  were  possible  because  the 
newer  version  of  the  CTH  code,  discussed  later,  was  more  stable  for  this  class  of  problems. 
Those  computations  are  fully  coupled  in  the  sense  that  the  detonation  of  the  munitions  stack, 
the  blast  loading  on  and  response  of  the  barricade;  and  loading  from  all  sources  on  and  the 
response  of  acceptor  stack  are  modeled  in  a  single,  continuous  computation.  The  first  study 
of  the  effects  of  standoff  variation  on  the  loading  and  response  of  the  acceptor  stack  was  done 
for  a  massive  water  barricade  having  a  trapezoidal  cross  section.^  The  barricade  is  identical 
to  that  modeled  in  the  uncoupled  pair  of  computations.^  That  first  study  also  showed  that 
the  sloping  sides  of  the  trapezoidal  water  barricade  were  effective  in  deflecting  air  blast  and 
explosive  products  upward  and  away  from  the  acceptor  stack.  The  normalized  blast  loading 
on  the  trapezoidal  barricade  and  its  whole-body  response  were  relatively  weak  functions  of 
the  inverse  of  the  normalized  standoff  distance.  The  blast  and  impact  loading  on  the  acceptor 
stack  was  a  three-stage  process.  The  the  first  stage  was  from  the  air  blast,  and  the  next  two 
stages  came  first  from  a  water  wave  at  the  top  of  the  barricade  and  second  from  the  impact 
of  the  lower  section  of  the  barricade.  The  blast  loading  on  the  acceptor  stack  was  negligible 
compared  to  that  from  the  impact  of  the  water.  The  normalized  impact  loading  of  the  water 
barricade  on  the  acceptor  stack  had  a  nearly  one-to-one  correspondence  with  the  inverse  of 
the  normalized  standoff  distance.  Peak  pressures  on  the  acceptor  stack  surface  facing  the 
barricade  were  below  300  MPa  (3  kbar),  considerably  less  than  the  500  MPa  (5  kbar)  in  the 
uncoupled  computations.  The  five  standoff  distances  in  the  first  fully  coupled  computational 
study^  were  3.048  m  (10.0  ft),  2.75  m  (9.02  ft),  2.50  m  (8.20  ft),  2.25  m  (7.38  ft),  and  2.00  m 
(6.56  ft). 

The  second  series  of  computations  modeling  the  effects  of  standoff  variation  on  the 
loading  and  response  of  the  acceptor  stack  was  done  for  a  relatively  thin  (1.17-m  width) 
water  barricade  having  a  rectangular  cross  section. ^  It  showed  that  the  thin  rectangular 
water  barricade  was  effective  in  deflecting  blast  upward  and  away  from  the  acceptor  stack,  but 
was  much  less  effective  in  keeping  explosive  products  from  impinging  on  the  acceptor  stack. 
The  computations  demonstrated  a  relatively  weak  inverse  functional  relationship  between 
normalized  values  of  the  standoff  distance  and  the  loading  on  and  whole-body  response 
of  the  barricade.  Similar  results  for  both  standoff  and  face  separation  were  found  for  the 
loading  on  and  whole-body  response  of  the  acceptor  stack,  except  for  a  stronger  functional 
relation  of  acceleration.  The  impact  loadings  on  the  acceptor  stack  bj^  the  thin  rectangular 
water  barricade  are  much  more  severe  than  those  reported^  for  the  massive  trapezoidal 
water  barricade.  Peak  pressures  on  the  acceptor  stack  were  approximately  2  GPa  (20  kbar), 
high  enough  to  represent  a  threat  of  initiating  a  chemical  reaction  in  munitions  within  the 
acceptor  stack. 

This  report  describes  a  series  of  fully  coupled  computations  for  the  same  three  standoff 
distances  as  were  evaluated  in  the  study  involving  the  thin  rectangular  water  barricade. 
Computation  980918  is  for  a  3.048-m  (10.0-ft)  standoff  (hereinafter  rounded  to  3.05  m  for 
simplicitv,  except  when  specifically  used  to  calculate  a  parameter).  Computation  9b0923  is 
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for  a  2.50-m  (8.20-ft)  standoff,  and  Computation  980924  is  for  a  2.00-m  (6.56-ft)  standoff. 
The  donor  and  acceptor  stacks  are  modeled  in  a  way  that  is  identical  to  that  used  in  the 
previous  fully  coupled  studies^'  ^  and  the  first  uncoupled  study.^  The  1.70-m-thick  barricade 
is  modeled  as  a  simple  rectangle  having  the  same  height  as  both  the  donor  and  acceptor 
stacks,  and  therefore  the  same  height  as  the  trapezoidal  and  thin  rectangular  barricades  in 
the  previous  studies. 


2.  COMPUTATIONAL  APPROACH  AND  GEOMETRY 

2.1.  The  Hydrocode  Model 

The  three  coupled  computations  that  are  reported  here  were  performed  using  the  then- 
latest  general-release  version,  CTH_9801,  of  the  CTH^  hydrocode  developed  at  Sandia  Na¬ 
tional  Laboratories  (SNL).  It  also  includes  the  May  1998  and  August  1998  “patches”  (i.e., 
coding  updates)  that  were  released  by  SNL.  CTH  solves  the  inviscid  Euler  equations  using 
a  second-order  accurate,  explicit  time-stepping  method.  A  brief  description  of  the  CTH 
hydrocode  was  given  in  a  previous  report.^  The  reader  is  referred  to  McGlaun  et  al.^  for  a 
full  discussion  of  the  CTH  hydrocode,  and  to  the  appropriate  users’  manuals  for  practical 
information  about  the  structure  and  use  of  the  CTHGEN^  grid  generation  code,  the  CTH® 
hydrocode,  and  their  supporting  utilities. 

The  three  computations  presented  here  were  performed  using  the  two-dimensional  (2-D) 
Cartesian  coordinates  system  option  in  CTH,  just  as  was  done  for  the  previous  computations. 
The  choice  of  2-D  Cartesian  coordinates  meant  that  the  computations  provided  a  worst-case 
blast  loading  for  the  simplified,  uncased  charge  of  condensed  high  explosives  by  eliminating 
the  possibility  of  having  any  compression  or  expansion  waves  in  the  direction  of  depth  of 
the  munitions  stacks  and  barricade.  Depth  is  a  measure  parallel  to  both  the  ground  and  the 
side  walls  of  the  munitions  stack,  and  normal  to  the  page  in  the  flow  field  plots  shown  later. 
In  effect,  the  donor  and  acceptor  stacks  and  the  barricade  have  an  infinite  depth  in  that 
coordinate  system.  In  the  CTH  hydrocode  model,  which  uses  the  centimeter-gram-second 
(cgs)  units  system,  this  implies  a  unit  depth  of  1.0  cm.  The  same  gridding  was  used  in 
all  computations  for  all  of  these  studies.  The  nominal  computational  cell  dimensions  are 
4.0  cm  in  both  Aa:  (width)  and  Ay  (height).  These  computations  were  performed  on  the 
Silicon  Graphics,  Inc.,  (SGI)  Origin  2000  unclassified  computers  at  the  ARL  Major  Shared 
Resource  Center  (MSRC)  at  Aberdeen  Proving  Ground  (APG),  MD.  This  is  one  of  four 
MSRCs  in  the  United  States  that  are  administered  by  the  High  Performance  Computing 
Modernization  Office  (HPCMO).  Each  of  these  2-D  Cartesian  computations  took  about 
250,000  central-processor-unit  (cpu)  seconds  (nominally,  70  hours).  Each  used  approximately 
150,000  flow  field  cells,  with  43  variables  per  cell  describing  the  materials  and  their  dynamic 
and  thermodynamic  attributes. 
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2.2.  The  Donor  Munitions  Stack 

The  donor  stack  is  modeled  in  the  same  way  as  in  the  uncoupled^  and  coupled^’  ^  studies 
described  in  previous  reports:  as  an  uncased  charge  with  no  packing  materials.  This  reduced 
the  analysis  to  one  of  blast  loading  only,  with  no  production  of  fragments  or  other  debris. 
The  explosive  mass  of  the  donor  stack  is  modeled  as  a  single,  condensed  charge  rather 
than  as  a  distributed  set  of  smaller  condensed  charges.  The  choice  of  the  munitions  in 
the  donor  munitions  stack  was  made  by  consulting  a  previous  ARL  report  on  fragment 
propagation  probabilities  by  Starkenberg  et  al.’^  The  donor  munitions  stack  was  assumed  to 
be  of  the  same  physical  dimensions  as  one  consisting  of  72  pallets  of  M107  155-mm  projectiles, 
stacked  three  pallets  high  by  four  wide  by  six  deep.  Each  pallet  contains  eight  rounds.  The 
dimensions  of  this  particular  stack  were  2.44  m  high  by  2.94  m  wide  by  2.19  m  deep  (8.00  ft 
by  9.63  ft  by  7.20  ft).  A  single  M107  round  can  contain  either  6.62  kg  (14.6  Ibm,  where 
“Ibm”  denotes  pounds  mass,  avoirdupois)  of  TNT  or  6.98  kg  (15.4  Ibm)  of  Composition-B 
(hereinafter  referred  to  as  “Comp-B”).  The  total  mass  of  a  pallet,  including  packaging,  is 
362  kg  (797  Ibm).^  Thus,  a  presumed  stack  of  72  pallets  of  M107  munitions  would  contain 
576  rounds,  having  a  total  mass  of  Comp-B  equal  to  4,024  kg  (8,870  Ibm).  For  simplicity, 
the  nominal  explosive  mass  of  Comp-B  for  this  computational  study  was  taken  as  4,000  kg 
(8,818  Ibm)  for  the  donor  stack,  the  regulatory  maximum.^  The  total  mass  of  an  actual 
stack  containing  72  pallets  of  M107  rounds  is  26,029  kg  (57,384  Ibm),  including  all  packaging 
materials.  This  equates  to  a  mass  of  118.61  kg/cm  of  depth  for  the  actual  stack  with  all 
materials  considered.  The  acceptor  stack  was  assumed  to  be  of  the  same  physical  dimensions 
and  total  mass  as  those  of  the  donor  stack. 

The  explosive  modeled  was  Comp-B,  taken  at  its  reference  density  of  1.72  g/cm^  in  its 
undetonated  state,  and  modeled^O  within  the  Sesamell  equation-of-state  (EOS)  package. 
The  SNL  Sesame  EOS  package  includes  tabular  data  for  high  explosives  and  separate  imple¬ 
mentations  of  data  for  the  Mie-Gruneisen,  Jones- Wilkins- Lee  (JWL),  and  ideal-gas  EOSs. 
The  explosive  charge  was  placed  within  the  computational  flow  field  with  its  center  coinci¬ 
dent  with  the  geometric  center  of  the  M107  donor  stack  described  before.  After  assigning 
the  donor  stack  the  nominal  explosive  mass  of  4,000  kg  and  using  the  actual  stack  depth  of 
2.19  m,  this  equated  to  an  explosive  charge  mass  of  approximately  18.227  kg/cm  of  depth 
of  the  stack  to  be  modeled  in  the  unit-depth  2-D  Cartesian  coordinates  flow  field  in  CTH. 
This  mass  of  Comp-B  was  modeled  as  a  rectangle  whose  width  and  height  are  in  direct  pro¬ 
portion  to  those  for  the  donor  stack.  Specifically,  the  explosive  charge  is  93.91  cm  high  and 
113.04  cm  wide.  This  is  the  full  width,  and  not  the  one-half  width  used  because  of  charge 
symmetry  about  the  left  boundary  in  the  2-D  computation.  The  charge  was  located  with 
its  center  of  mass  121.92  cm  above  the  ground  plane.  The  ground  plane  was  designated  as 
a  frictionless,  perfectly  reflective  boundary. 

A  small  central  section  of  the  explosive  charge  at  the  left  symmetry  boundary  served 
as  a  computational  “booster”  charge.  It  was  detonated  using  the  programmed  burn'^  model 
using  a  constant  detonation  velocity  7.98  km/s  for  reference-density  Cornp-B.^^  This  model 
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simulates  the  complete  detonation' of  any  part  of  an  explosive  charge  that  is  passed  by  the 
expanding  theoretical  detonation  front  moving  at  that  constant  velocity.  The  remainder  of 
the  detonation  was  modeled  using  the  “history  variable  reaction  burn”  (HVRB)  model.^l 
The  HVRB  model  evaluates  the  thermodynamic  state  of  a  mass  of  undetonated  explosive 
in  a  given  computational  flow  field  cell  to  determine  if  that  material  should  be  numerically 
“reacted”  to  simulate  its  detonation  in  that  time  step.  The  detonation  initiation  point  was 
located  at  the  center  of  the  explosive  charge  at  the  (X,Y)  point  (0.0,  121.92  cm)  on  the  left 
symmetry  boundary. 


2.3.  The  Barricade 

The  barricade  shape  chosen  for  the  computations  was  a  simple  rectangle  with  a  nominal 
height  of  2.44  m  (8.0  ft).  This  is  the  same  height  as  that  for  the  munitions  stacks  and  for  the 
trapezoidal  and  thin  rectangular  barricades  studied  previously.^’  ^  The  width  chosen  for  this 
thick  rectangular  barricade,  1.70  m,  is  equal  to  the  width  of  the  thin  rectangular  barricade, 
1.17  m,  times  a  factor  of  1.45.  The  factor  of  1.45  was  chosen  because  it  was  numerical  value 
of  the  density  of  the  sand,  1.45  g/cm^,  in  a  sand-filled  Concertainer  barricade  that  was  tested 
against  a  detonated  pallet  of  M107  rounds  for  ARL^^-  by  the  Aberdeen  Test  Center  (ATC). 
(Concertainer  is  marketed  by  Hesco  Bastion  Limited. )  This  produced  a  water  barricade 
with  the  same  total  mass  as  one  made  out  of  sand  with  a  density  of  1.45  g/cm^  and  a  1.17- 
m  width.  This  was  done  because  earlier  attempts  to  model  a  rectangular  sand  barricade 
matching  the  Concertainer  barricade  failed  because  of  numerical  stability  problems  soon 
after  the  initial  interaction  of  the  blast  wave  with  the  barricade.  These  computations  with 
the  thick  rectangular  water  barricade  were  performed  as  a  compromise  to  at  least  simulate 
the  same  mass  and  inertial  effects  as  a  sand-filled  Concertainer  barricade  with  a  height  of 
2.44  m.  Other  than  using  this  nominal  value  of  1.70  m  for  the  width  of  the  rectangular  water 
barricade,  the  rectangular  water  barricade  simulated  in  the  computations  reported  here  has 
no  relation,  direct  or  indirect,  to  Concertainer  or  any  other  product  by  Hesco  Bastion.  Any 
data  related  to  or  evaluation  of  the  rectangular  water  barricade  simulated  here  also  should 
not  be  construed  as  having  any  relation  to  any  commercial  product  by  Hesco  Bastion.  An 
idealized  cross  section  that  has  no  internal  air  spaces  and  consists  only  of  water  is  assumed. 
No  construction  or  supporting  materials  are  considered.  The  mass  of  water  for  the  barricade 
is  41.49  kg/cm  of  depth.  The  water  in  the  barricade  was  modeled  using  the  CTH  Sesame 
EOS  for  water.  The  face  of  the  barricade  closest  to  the  donor  stack  was  placed  at  the 
defined  standoff  distance,  which  was  varied  from  3.05  m  (10  ft)  to  2.00  m  (6.56  ft),  from  the 
nearest  side  of  the  donor  stack.  The  standoff  distance  in  the  computations  here  is  measured 
from  the  face  of  what  would  have  been  the  actual  side  of  the  donor  munitions  stack,  not  the 
condensed  explosive  charge  representing  the  stack. 
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2.4.  The  Acceptor  Munitions  Stack 

The  acceptor  munitions  stack  was  modeled  in  all  computations  as  a  simple,  relatively 
inert  mass  of  ironl”^  with  the  same  height  (2.44  m)  and  width  (2.94  m)  as  the  reference  M107 
munitions  stack.  This  is  identical  to  the  way  the  acceptor  stack  was  modeled  in  all  previous 
uncoupled^  and  coupled^’  ^  computations.  The  acceptor  stack  in  each  computation  was  lo¬ 
cated  at  a  standoff  distance  equal  to  that  between  the  donor  stack  and  the  barricade.  The 
purpose  in  modeling  the  acceptor  stack  as  a  full-sized  mass  of  iron  was  for  the  convenience  of 
having  a  massive,  relatively  non-responding  object  with  the  correct  physical  dimensions  in 
order  to  observe  wave  interactions  on  the  surface  and  to  provide  surface  blast  loading  data 
through  the  use  of  CTH’s  “tracer”  particles  placed  in  the  air  near  the  surfaces.  Tracer  parti¬ 
cles  are  massless  points  that  are  specified  at  desired  locations  by  the  user  at  grid  generation 
time.  They  may  be  fixed  in  computational  space  or  be  free  to  move  along  one  or  more  of 
the  principal  axes  in  the  grid.  A  relatively  full  complement  of  data  describing  the  thermody¬ 
namic  state  and  other  physical  parameters  at  the  location  of  each  tracer  is  recorded  for  later 
processing  by  the  user.  W  hen  analyzing  the  whole-body  response  of  the  acceptor  stack  later 
in  this  report,  the  correct  acceptor  stack  mass  (118.61  kg/cm  of  depth)  was  used  to  compute 
the  motion  of  the  acceptor  stack  from  the  X-direction  momentum  of  the  massive  iron  stack. 
Some  of  the  details  of  the  blast  development  in  this  computational  series  and  subsequent 
interactions  between  the  blast  and  the  barricade  and  then  the  barricade  and  the  acceptor 
stack  are  surely  artifices  of  the  simplified  geometries,  but  the  overall  dynamics  appear  to  be 
quite  reasonable. 


3.  THE  HYDROCODE  COMPUTATIONS 
3.1.  Flow  Field  Development 

Computation  980918  simulated  a  fully  coupled  blast  and  impact  loading  sequence  at  a 
standoff  of  3.05  m  (10.0  ft).  Even  though  the  standoff  is  the  same  at  that  in  Computation 
980505^  and  Computation  980825,'^  there  are  significant  differences  in  the  barricades  and 
their  relative  positions  beyond  the  obvious  differences  in  the  cross  sections  and  slopes  of  the 
left  and  right  faces.  The  trapezoidal  water  barricade  mass  is  58.71  kg/cm  of  depth  with  an 
X-direction  distance  of  its  center  of  mass  equal  to  4.96  m  from  the  right  face  of  the  donor 
stack.  The  thin  rectangular  water  barricade  has  a  mass  of  28.61  kg/cm  of  depth  with  an 
X-direction  distance  of  its  center  of  mass  equal  to  3.63  m  from  the  right  face  of  the  donor 
stack.  The  thick  rectangular  water  barricade  has  a  mass  of  41.49  kg/cm  of  depth  with  an 
X-direction  distance  of  its  center  of  mass  equal  to  3.90  m  from  the  right  face  of  the  donor 
stack.  These  parameters  are  summarized  in  Table  1  in  dimensional  form  and  also  normalized 
by  dividing  by  the  particular  value  for  the  massive  trapezoidal  water  barricade  for  the  3.05-m 
standoff  as  well  as  the  other  standoff  distances. 
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Table  1.  Barricade  Geometries  and  Comparisons  for  Three  Standoffs. 


Geometry 

Massive 

Trapezoidal 

Thick 

Rectangular 

Thin 

Rectangular 

Mass 

(kg/ cm  of  depth) 

58.71 

41.49 

28.61 

Normalized  Mass 
(-) 

1.000 

0.7067 

0.4873 

Inside  Angle  of 
Side  to  the 
Vertical 
(Degrees) 

30.0 

0.0 

0.0 

3.048-m  Standoff, 
Donor  Right  Face 
to  Barricade 
Center  of  Mass 
(m) 

4.9558 

3.8988 

3.6347 

Normalized  by 
Massive  Trapezoidal 
Value 
(-) 

1.0000 

0.7867  . 

0.7334 

2.50-m  Standoff, 
Donor  Right  Face 
to  Barricade 
Center  of  Mass 
(m) 

4.4078 

3.3508 

3.0867 

Normalized  by 
Massive  Trapezoidal 
Value 
(-) 

1.0000 

0.7602 

0.7003 

2.00-m  Standoff, 
Donor  Right  Face 
to  Barricade 
Center  of  Mass 
(m) 

3.9078 

2.8508 

2.5867 

Normalized  by 
Massive  Trapezoidal 
Value 
(-) 

1.0000 

0.7295 

0.6619 
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Figure  1  shows  the  computational  flow  field  at  the  start  of  Computation  980918.  This 
is  the  instant  of  the  initiation  of  the  detonation  (hereinafter  referred  to  as  “initiation”), 
with  time  defined  to  be  equal  to  zero.  The  “Y”  axis  at  the  left  of  the  figure  represents  the 
height  measured  from  the  ground  plane.  In  this  simple  2-D  Cartesian  coordinate  system, 
the  left  boundary  at  the  Y  axis  is  designated  as  a  frictionless,  perfectly  reflective  plane  of 
symmetry.  The  “X”  axis  represents  the  measure  of  width  in  the  system  and  coincides  with 
the  frictionless,  perfectly  reflective  ground  plane.  The  X  direction  therefore  represents  the 
measure  of  thickness  of  the  barricade.  The  Y  axis  at  the  X  =  0.0  location  is  also  a  vertical 
bisector  of  the  donor  stack.  The  air  in  the  flow  field,  modeled  with  data  from  Graboske^^ 
within  the  Sesame^^  EOS,  is  shown  with  the  color  yellow.  The  top  and  right  transmissive 
boundaries  are  marked  by  the  top  and  right  edges  of  that  yellow  region.  These  transmissive 
boundaries  were  designated  as  zero-gradient,  outflow-only  boundaries  to  minimize  the  pos¬ 
sibilities  of  generating  spurious,  mathematically  generated  reflected  waves  or  inflows  when 
those  boundaries  are  struck  by  large-gradient  outflows.  The  explosive  charge  representing 
the  donor  stack  is  shown  as  the  red  (one-half)  rectangle  on  the  left  symmetry  boundary,  the 
water  barricade  is  shown  as  the  blue  rectangle,  and  the  acceptor  stack  is  shown  as  the  black 
rectangle  (the  object  closest  to  the  right  transmissive  boundary).  These  settings  and  general 
descriptions,  except  for  the  shape  of  the  barricade,  are  the  same  as  those  for  the  previous 
series^  2  of  computations.  In  order  to  facilitate  direct  comparison  of  the  flow  fields  for  the 
3.05-m  standoff  for  the  thick  rectangular  versus  the  thin  rectangular  and  trapezoidal  barri¬ 
cades,  representative  flow  fields  from  Computations  980505^  and  980825^  are  also  presented. 
Figure  2  shows  the  flow  field  at  the  instant  of  initiation  for  Computation  980505  for  the 
trapezoidal  barricade.  Figure  3  shows  the  flow  field  at  the  same  initial  time  for  Computa¬ 
tion  980825.  A  comparison  between  Figures  1,  2,  and  3  provides  a  good  visual  indication  of 
the  differences  in  spacing  and  mass  of  the  different  barricades  for  the  same  standoff  distance. 
The  spatial  extent  of  the  flow  field  in  each  direction  is  the  same  for  all  three  computations, 
as  is  the  fineness  of  the  computational  gridding.  The  donor  stack,  its  initiation  point,  and  its 
location  are  also  identical.  The  acceptor  stack  itself  is  the  same  in  all  three  computations, 
but  its  X  location  varies  because  of  the  differences  in  the  width  of  the  barricades  at  their 
bases. 

Figure  4  shows  the  computational  flow  field  for  980918  at  5.0  ms  after  the  initiation  of  the 
donor  stack.  The  detonation  process  had  already  been  completed  by  this  time  (theoretically 
at  0.092  ms).  The  expanding  explosive  products  and  leading  shock  have  deformed  and 
accelerated  the  barricade.  A  section  of  the  expanding  explosive  products  has  already  passed 
over  the  acceptor  stack.  The  lower  section  of  the  barricade  is  translating  laterally  toward  the 
acceptor  stack  as  a  relatively  unified  block  with  the  lower  section  of  the  barricade  leading  the 
rest  of  the  barricade.  The  bottom-most  part  of  the  barricade  is  approximately  1.4  m  away 
from  the  acceptor  stack  left  face.  There  is  some  shearing  off  of  the  top  of  the  barricade.  As 
yet,  no  part  of  the  barricade  has  arrived  at  the  left  face  of  the  acceptor  stack.  No  significant 
amounts  of  explosive  products  appear  to  have  reached  the  acceptor  stack.  For  comparison. 
Figure  5  shows  the  computational  flow  field  for  980505  at  5.0  ms  after  the  initiation  of  the 
donor  stack.  The  more  massive  trapezoidal  water  barricade  in  that  computation  is  also 


8 


2DR  Block  1  X  (m) 

980505  TRAP  H20  BARR,  3.05  M  STANDOFF 

FPQDPl  G  06/16/98  16:40:26  CTHGEN  0  Time=0. 


figure  2.  Flow  Field  at  Time  =  0.0  for  Computation  980505,  ;3.05-m  StaridofF,  Massive 
Frapezoidal  Barricade. 
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Figure  ;3.  Flow  Field  at  Time  =  0.0  for  Computation  980825,  3.1 
angular  Barricade. 
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showing  the  leading  action  of  the  lower  section  of  the  barricade,  but  with  significantly  less 
motion.  This  leading  action  is  occurring  in  a  two-fold  manner.  First,  the  base  of  the  barricade 
is  simply  translating  laterally  toward  the  acceptor  stack  in  response  to  the  blast  loading. 
Second,  "there  is  a  single  wave  developing  on  the  lower-right  surface  of  the  barricade  that  is 
already  beginning  to  lead  the  rest  of  the  lower  section  of  the  barricade  and  travel  up  the  right 
surface  of  the  barricade.  At  this  time,  the  tip  of  that  wave  is  approximately  2.5  m  from  the 
left  surface  of  the  acceptor  stack  and  0.4  m  above  the  ground  plane.  This  wave  development 
was  discussed.in  detail  previously.^  Conversely,  Figure  6  shows  the  computational  flow  field 
for  980825  for  the  thin  rectangular  barricade  at  5.0  ms.  It  shows  the  greatest  distortion  and 
movement  toward  the  acceptor  stack.  The  base  of  the  barricade  is  approximately  1.1  m  fiom 
the  left  surface  of  the  acceptor  stack.  It  is  leading  the  rest  of  the  barricade,  except  for  the 
shearing  and  dispersion  at  its  top. 

Figure  7  shows  the  computational  flow  field  for  the  thick  rectangular  barricade  at  10.0  ms 
after  initiation.  The  impact  of  the  barricade  on  the  left  surface  of  the  acceptor  stack  is  nearly 
at  the  end  of  its  full-interaction  phase,  covering  the  entire  left  surface  of  the  acceptor  stack. 
The  lower  part  of  the  barricade  is  still  in  the  form  of  a  relatively  thick  wedge  of  water  on 
the  acceptor  stack  left  surface.  The  upper  part  of  the  barricade  has  distorted  an  additional 
amount  and  is  extending  farther  upward.  Explosive  products  are  approaching  the  top-rear 
corner  of  the  acceptor  stack.  The  interactions  for  the  other  two  computations  at  the  same 
time  are  in  very  different  phases  from  this  and  from  one  another.  Figure  8  for  980505  shows 
that  the  upward-moving  wave  on  the  right  surface  of  the  barricade  is  continuing  to  develop. 
The  base  of  the  barricade  is  approximately  2.0  m  away  from  the  left  surface  of  the  acceptor 
stack  The  tip  of  the  wave  is  now  1.2  m  above  the  ground  plane  and  1.2  rn  away  from  the 
acceptor  stack  left  surface.  Figure  9  for  the  thin  rectangular  barricade  in  980825  has  already 
completed  its  impact  on  the  acceptor  stack,  has  rebounded  Irom  it,  and  is  moving  in  the 
negative  X  direction.  Although  it  still  has  appearance  of  a  continuous  structure  in  front  of 
the  acceptor  stack  left  surface,  it  is  quite  ragged.  The  remainder  of  the  banncade  has  been 
dispersed  upward  and  away  from  the  acceptor  stack.  Some  explosive  products  are  relatively 
close  to  the  top-rear  corner  of  the  acceptor  stack. 

Figure  10  shows  the  flow  field  for  the  thick  rectangular  barricade  at  15.0  ms.  By  this 
time,  the  barricade  has  rebounded  from  the  acceptor  stack  left  face,  but  still  shows  a  mod¬ 
erate  amount  of  coherent  structure  in  front  of  the  acceptor  stack.^  Much  of  the  barricade  has 
been  dispersed  upward  to  two  or  more  multiples  of  its  original  height  as  well  as  being  spread 
downstream.  Explosive  products  are  approaching  the  top  face  of  the  acceptor  stack.  Fig¬ 
ure  11  shows  the  flow  field  at  15.0  ms  for  the  trapezoidal  barricade  in  Computation  980o0d. 
It  shows  a  distorted  but  still-intact  barricade  with  the  tip  of  the  wave  on  its  right-rear  face 
almost  at  the  point  of  its  first  contact  with  the  left  face  of  the  acceptor  stack.  The  baae  of 
the  barricade  is  1.4  m  away  from  the  left  surface  of  the  acceptor  stack.  Most  of  the  air  blast 
and  virtually  all  of  the  explosive  products  have  been  deflected  upward  and  away  from  the 
acceptor  stack  by  the  trapezoidal  barricade.  Figure  12  shows  that  by  15.0  ms  after  initiation, 
the  thin  rectangular  barricade  no  longer  has  any  useful  structural  integrity  and  no  readily 
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Figure  4.  Flow  Field  at  Time  =  5.0  ms  for  Computation  980918,  3.05-m  Standoff,  Thick 
Rect  angular  Barri  cade. 
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Figure  5.  Flow  Field  at  Time  =  5.0  ms  for  Computation  980505.  3.05-m  Standoff,  Massive 
Trapezoidal  Barricade. 
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Figure  6.  Flow  Field  at  Time  =  5.0  ms  for  Computation  980825,  3.05-m  Standoff,  Thin 
Rect  angular  B  arricade. 


2DR  Block  1  X 

980918  RECT  (1.45 w)  H20  BARR,  3.05  M  STANDOFF 
IRPELA  09/19/98  21:49:13  CTH  5575  Time=1.00016x10- 


Figiire  7.  Flow  Field  at  Time  =  10.0  ms 
Rectangular  Barricade. 
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Figure  8.  Flow  Field  at  Time  =  10.0  ms  for  Computation  980505,  3.05-m  Standoff,  Massive 
Trapezoidal  Barricade. 
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Figure  9.  Flow  Field  at  Time  =  10.0  ms 
Rectangular  Barricade. 


for  Computation  980825,  3.05-m  Standoff,  .riiin 


recognized  shape.  A  small  amount  of  explosive  products  appears  to  be  making  contact  with 
the  lower-left  corner  of  the  acceptor  stack,  and  a  large  region  of  explosive  products  seems  to 
be  moving  very  close  to  the  top-left  corner  of  the  stack. 

Figure  13  shows  the  computational  flow  field  for  the  thick  rectangular  barricade  at 
20.0  ms  after  initiation.  The  rebound  of  the  barricade  from  the  left  surface  of  the  acceptor 
stack  is  continuing,  and  the  barricade  is  in  the  process  of  being  dispersed  into  the  rest  of 
the  flow  field.  Explosive  products  seem  to  now  be  contacting  at  least  part  of  the  top  suiface 
of  the  acceptor  stack  and  are  close  to  the  bottom-left  corner.  Figure  14  for  980505  shows 
the  first  interaction  of  water  from  the  wave  on  the  right  face  of  the  trapezoidal  barricade 
with  the  top  section  of  the  acceptor  stack.  The  bottom  of  the  barricade  is  approximately 
0.9  m  from  the  acceptor  stack  left  surface.  The  barricade  is  still  largely  intact,  providing 
good  protection  for  the  acceptor  stack,  and  no  explosive  products  are  near  the  surfaces  of 
the  acceptor  stack.  Figure  15  shows  the  computational  flow  field  at  20.0  ms  after  initiation. 
Much  of  the  top  face  of  the  barricade  has  contact  or  near  contact  with  explosive  products, 
as  does  its  lower-left  corner.  Some  explosive  products  are  nearing  the  back  face.  There  is 
no  longer  any  structurally  meaningful  section  of  barricade  in  front  of  the  left  surface  of  the 
acceptor  stack  or  anywhere  else  in  the  flow  field. 

Progressing  further  in  time.  Figures  16  and  17  show  the  computational  flow  field  for 
the  thick  rectangular  barricade  at  30.0  ms  and  40.0  ms,  respectively,  after  initiation.  They 
show  remnants  of  the  barricade  being  pressed  once  again  against  the  left  surface  of  the 
acceptor  stack  and  then  dispersed  into  the  flow  field.  By  40.0  ms,  there  is  general  contact  of 
explosive  products  with  the  acceptor  stack  surfaces.  There  are  no  recognizable,  contiguous 
sections  of  the  barricade  left  in  the  computational  flow  field.  A  significant  portion  of  the 
original  barricade  mass  has  exited  the  flow  field.  In  contrast,  the  interaction  of  the  massive 
trapezoidal  barricade  with  the  left  surface  of  the  acceptor  stack  had  not  yet  reached  its  peak 
by  30  ms.  The  base  was  still  about  0.5  m  from  the  left  surface  of  the  acceptor  stack.  By 
40  ms,  the  interaction  had  passed  its  peak  and  the  lower  section  of  the  barricade  was  in  a 
rebound  phase  and  showing  a  relatively  good  continuous,  though  highly  distorted,  structure 
in  front  of  the  acceptor  stack.  Those  flow  fields  are  not  shown  here,  but  may  be  seen  in  the 
report  on  the  trapezoidal  barricade  computations.^  The  flow  fields  for  the  thin  rectangular 
barricade  at  30  and  40  ms  show  that  the  barricade  has  been  effectively  removed  from  the 
flow  field  and  that  explosi  ve  products  are  in  general  contact  with  the  acceptor  stack.  A  more 
complete  description  along  with  plots  of  those  flow  fields  may  be  seen  in  the  report  on  the 
thin  rectangular  barricade  computations.'^ 

From  a  qualitative  point  of  view  from  comparing  the  flow  fields  for  the  3.05-m  standoff, 
there  is  a  clear  hierarchy  in  the  effectiveness  of  these  three  different  water  barricades.  The 
thick  rectangular  barricade  was  more  effective  for  a  longer  period  of  time  than  was  the  thin 
rectangular  barricade.  However,  both  of  the  rectangular  barricades  were  considerably  less 
effective  in  protecting  the  acceptor  stack  than  was  the  more  massive  trapezoidal  barricade. 
This  statement  is  quantified  in  the  following  sections  of  this  report.  First,  sets  of  similai 
sequences  of  the  flow  fields  for  the  thick  rectangular  water  barricade  for  a  2.50-m  and  a 
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Figure  10.  Flow  Field  at  Time  =  15.0  ms  for  Computation  980918,  3.05-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  11.  Flow  Field  at  Time  =  15.0  ms  for  Computation  980505,  3.05-m  Standoff,  Massive 
Trapezoidal  Barricade. 
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Figure  12.  Flow  Field  at  Time  =  15.0  ms  for  Computation  980825,  3.05-Tn  Standoff,  Thin 
Rectangular  Barricade. 
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Figure  13.  Flow  Field  at  Time  =  20.0  ms  for  Computation  980918,  3.05-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  14.  Flow  Field  at  Time  =  20.0  ms  for  Computation  980505,  3.05-m  Standoff,  Massive 
Trapezoidal  Barricade. 
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Figure  15.  Flow  Field  at  Time  =  20.0  ms  for  Computation  980825,  3.05-m  Standoff,  Tlnn 
Rectangular  Barricade. 
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Figure  16.  Flow  Field  at  Time  =  30.0  ms 
Rectangular  Barricade. 


for  Computation  980918,  3.05-m  Standoff,  J  hiclc 


2  00-m  standoff  are  shown.  They  both  show  similar  behavior  to  that  for  980918  so  fewer 
snapshots  in  time  are  presented,  and  no  direct  comparison  to  corresponding  flow  fields  tor 
the  previous  two  series  of  computations!-  2  the  massive  trapezoidal  and  thin  rectangular 

barricades,  respectively,  are  shown. 


Computation  980923  simulated  a  standoff  distance  of  2.50  m  for  the  same  thick  rect¬ 
angular  water  barricade  as  in  980918.  Figure  18  shows  the  computational  flow  field  at 
time  =  0.0.  Only  the  standoff  distance  has  been  changed  from  that  for  980918,  theieby 
moving  the  barricade  leftward  in  the  flow  field  toward  the  donor  stack  by  0.55  m  and  the 
acceptor  stack  leftward  by  1.1  m.  The  gridding  for  the  computational  flow  field,  the  spatial 
dimensions  of  the  flow  field,  the  munitions  stacks,  and  the  barricade  are  all  identical  to  those 
for  980918  (see  Figure  1).  Figure  19  shows  the  computational  flow  field  at  time  =  5.0  ins 
for  Computation  980923.  This  shows  similar  behavior  to  that  shown  in  Figure  4  for  980918, 
except  that  the  barricade  is  closer  in  space  and  time  to  its  initial  contact  with  the  acceptor 
stack  left  face.  The  base  of  the  barricade  is  0.7  m  from  the  left  surface  of  the  acceptor  stack 
and  leading  the  main  body  of  the  translating  barricade.  Explosive  products  are  approaching 
the  top-left  corner  of  the  acceptor  stack,  possibly  being  drawn  in  by  the  low-pressure  center 
of  a  vortex  originating  from  that  corner.  Figure  20  shows  the  flow  field  at  10.0  ms  after 
initiation.  This  is  approximately  2.0  ms  after  the  peak  pressures  have  been  incurred  on  the 
acceptor  stack  left  face.  This  is  documented  and  quantified  in  a  later  section  of  this  report. 
The  rebound  from  the  acceptor  stack  left  surface  of  the  water  from  the  bai-ricade  is  now 
under  way  The  section  of  the  barricade  in  front  of  the  acceptor  stack  still  has  a  reasonable 
structural  integrity.  Explosive  products  are  near  the  top-left  corner  of  the  acceptor  stack_ 
Figure  21  shows  the  flow  field  for  980923  at  15.0  ms.  The  rebound  from  the  left  surface  of 
the  acceptor  stack  of  the  water  is  continuing,  and  the  structural  integrity  of  that  section 
of  water  is  showing  rapid  degradation.  Figure  22  shows  that  by  20.0  ms,  any  protection  of 
the  acceptor  stack  by  the  barricade  has  essentially  ended.  Figures  23  and  24  for  30  ms  and 
40  ms,  respectively,  show  the  re.nina.nts  of  the  barricade  being  swept  out  of  the  flow  field  and 
explosWe  products  in  general  contact  with  the  acceptor  stack. 


The  last  computation  in  this  series  for  thick  rectangular  water  barricades,  Computation 
9809^4  simulated  a  standoff  distance  of  2.00  m.  As  was  the  case  for  Computation  980923, 
the  initial  flow  field  and  layout  of  objects  are  the  same  as  for  Computation  980918  except 
for  the  standoff  distance  and  the  subseciuent  leftward  shifting  of  the  barricade  and  acceptor 
stack  Figure  25  shows  the  computational  flow  field  at  time  =  0.0.  Figure  26  shows  t  le 
flow  field  for  Computation  980924  at  5.0  ms.  At  this  time  the  lower,  leading  section  of  the 
barricade  is  just  about  to  strike  the  left  face  of  the  acceptor  stack.  Explosive  products  are 
approaching  the  top-left  corner  of  the  acceptor  stack.  They  will  soon  be  pushed  upward 
and  awav  from  it  by  the  air  being  forced  upward  along  the  left  surface  of  the  st^^k  by  the 
water  impact  moving  up  the  stack.  Figure  27  shows  the  flow  field  for  980924  at  10-0  ms. 
The  rebound  from  the  acceptor  stack  left  surface  of  the  water  is  well  under  way.  There 
are  trace  amounts  of  explosive  products  near  the  top  surface  of  the  acceptor  stack,  just 
downstream  from  the  top-left  corner.  Figures  28,  29,  30,  and  31  show  subsequent  images  o 
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Figure  18.  Flow  Field  at  Time  =  0.0  for  Computation  980923,  2.30-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  19.  Flow  Field  at  Time  =  5.0  ms  for  Computation  980923,  2.30-m  Standoff, 
Rectangular  Barricade. 
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Figure  20.  Flow  Field  at  Time  =  10.0  ms  for  Computation  980923,  2.50-m  Standoff,  Thjck 
Rect angular  Barricade. 
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Figure  21.  Flow  Field  at  Time  =  15.0  ms  for  Computation  98092-3,  2.50-m  StaiidofF,  Thick 
Rectangular  Barricade. 
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Figure  22.  Flow  Field  at  Time  =  20.0  ms  for  Computation  980923,  2.50-m  Standoff,  Tliick 
Rectangular  Bari'icade. 
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Figure  23.  Flow  Field  at  Time  =  30.0  ms  for  Computation  980923,  2.50-m  Standoff,  T'liick 
Rectangular  Barricade. 


the  computational  flow  field  at  time  =  15.0  ms,  20.0  ms,  30.0  ms,  and  40.0  ms,  respectively, 
for  Computation  980924.  Collectively,  the  figures  show  the  same  qualitative  behavior  as 
those  for  Computations  980918  and  980923,  with  events  happening  at  a  faster  rate  with 
respect  to  time  simply  because  of  the  closer  initial  standoff. 
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Figure  '25.  Flow  Field  at  Time  =  0.0  for  Computation  980924,  2.00-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  26.  Flow  Field  at  Time  =  5.0  ms  for  Computation  980924,  2.00-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  27.  Flow  Field  at  Time  =  10.0  ms  for  Computation  980924,  2.00-m  Standoff,  Thick 
Eectangulax  Barricade. 
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Figure  28.  Flow  Field  at  Time  =  15.0  ms  for  Computation  980924,  2.00-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  29.  Flow  Field  at  Time  =  20.0  ms  for  Computation  980924,  2.00-m  Standoff,  Tliick 
Rectangular  Barricade. 
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Figure  30.  Flow  Field  at  Time  =  30.0  ms  for  Computation  980924,  2.00-m  Standoff,  Thick 
Rectangular  Barricade. 
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Figure  31.  Flow  Field  at  Time  =  40.0  ms  for  Computation  980924,  2.00-m  StaxidofT,  Thick 
Rectangular  Barricade. 


3.2.  Barricade  Dynamics 


Figure  32  shows  the  bulk  momentum  per  centimeter  depth  of  the  water  barricade  in 
the  X  direction.  Positive  momentuin  is  in  the  direction  of  increasing  values  of  X.  Figure  32 
includes  the  X-direction  momentum  for  each  of  the  three  fully  coupled  computations  in  this 
series  for  the  thick  rectangular  barricade,  plus  that  for  Computation  980b25  for  the  thin 
rectangular  water  barricade’^  and  Computation  980505  for  the  massive  trapezoidal  water 
barricadel  a  3.05-m  standoff.  Hereinafter,  any  use  of  the  term  “momentum”  or  the  other 
variables  (e.g.,  velocity,  acceleration,  and  displacement)  derived  from  it  should  be  construed 
as  referring  to  the  bulk  value  in  the  X  direction  per  centimeter  of  depth,  unless  specifically 
stated  otherwise.  The  term  “bulk”  is  implied  but  used  only  sparingly  in  order  to  avoid 
repetition.  The  momentum  shown  here  is  the  combined  momentum  for  all  of  the  water  in 
the  flow  field  at  each  computational  time  step.  Values  for  the  mass  and  momentum  for 
the  water  (and  all  other  materials)  are  saved  after  each  time  step.  During  any  given  time 
step  later  in  the  computations,  some  water  flows  out  of  the  flow  field  through  either  or 
both  of  the  top  and  right  transmissive  boundaries.  Each  of  the  three  curves  for  the  tliick 
rectangular  barricade  shows  a  very  rapid,  monotonic  initial  increase  in  momentum  with 
decreasing  standoff.  The  curve  for  Computation  980924  (2.00-m  standoff)  shows  the  first, 
very  abrupt  decrease  in  momentum  after  its  peak  at  5.1  ms,  followed  by  980923  (for  the 
2.50-m  standoff)  after  its  peak  at  6.6  ms,  and  then  by  980918  (3.05-m  standoff)  after  its 
peak  at  8.3  ms.  This  very  rapid  drop  in  each  curve  for  the  thick  rectangular  barricade  is 
because  of  the  nearly  simultaneous  bottom-to-top  strike  of  the  barricade  on  the  left  face  of 
the  acceptor  stack.  As  may  be  seen  in  Figure  32,  the  thick  rectangular  barricade  delivers 
most  of  its  momentum  to  the  acceptor  stack  over  a  very  short  period  of  time,  as  does  the 
representative  plot  for  the  thin  rectangular  barricade^  at  a  3.05-m  standoff,  labeled  “Thin 
Rect,  Standoff  3.05  m.”  The  fifth  curve  shown  in  Figure  32,  labeled  “Trap  Standoff  3.05^m,’ 
is  for  the  momentum  of  the  massive  trapezoidal  water  barricade  at  a  3.05-m  standoff.  It 
shows  a  more  gradual  increase  to  a  considerably  lower  peak  momentum  with  a  two-stage, 
much  smaller  total  decrease  in  momentum. 

It  is  interesting  to  point  out  here  that,  while  there  is  not  much  difference  in  both  the  peak 
momenta  and  the  values  to  which  the  momenta  for  the  both  the  thick  and  thin  rectangular 
barricades  decrease  by  10  ms  as  a  function  of  standoff  distance,  there  is  a  great  difference 
in  both  peak  momentum  and  change  in  momentum  that  can  be  seen  when  comparing  both 
rectangular  barricades  with  the  massive  trapezoidal  barricade  at  the  same  3.05-m  standoff 
distance.  Table  2  contains  a  summary  of  several  X-directio.n  parameters  that  describe  some  of 
the  bulk  motion  of  the  barricade  for  the  various  computations.  First  among  those  parameters, 
after  the  computation  numbers,  standoff  distances,  and  then  barricade  masses,  are  the  peak 
X-direction  bulk  momentum  values  for  the  barricade,  along  with  their  respective  times  of 
occurrence,  listed  with  more  significant  figures  than  were  tjTiically  used  in  the  text  for 
completeness.  In  order  to  facilitate  comparisons,  the  first  column  of  numbers  is  for  the 
massive  trapezoidal  barricade  at  a  3.05-m  standoff,  followed  by  data  for  the  thin  rectangular 
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barricade  at  a  3.05-m  standoff,  and  then  the  three  standoffs  for  the  thick  rectangular  bar¬ 
ricade.  The  rest  of  the  parameters  in  the  table  are  discussed  in  the  following  paragraphs. 


Figure  32.  Water  Barricade  X-Direction  Momentum  Toward  the  Acceptor  Stack,  Compu¬ 
tations  980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and 
980505  ( Trapezoi  dal) . 


After  each  computational  time  step,  both  the  total  momentum  and  mass  of  the  water 
in  the  flow  field  are  known.  The  X-direction  bulk  velocity  (hereinafter  referred  to  as  “X- 
direction  velocity”)  of  the  barricade  may  be  computed  for  eadi  time  step  by  diviclmg  the 
instantaneous  momentum  by  the  corresponding  mass.  The  X-direction  velocity  of  the  water 
barricade  toward  the  acceptor  stack  for  each  standoff  for  the  thick  rectangular  barricade, 
plus  the  velocity  for  the  trapezoidal  barricade  and  the  thin  rectangular  barricade  at  a  3.05-m 
standoff,  is  shown  in  Figure  33.  These  curves  are  essentially  scaled  variants  of  the  momentum 
curves  shown  in  Figure  32  and  therefore  show  the  same  relative  behavior  described  for 
the  momenta.  The  peak  X-direction  velocity  for  the  thick  rectangular  barricade  at  each 
standoff  distance  is  355.3  m/s  at  8.33  ms  (3.05-m  standoff),  355.9  m/s  at  6.56  ms  (2.50-m 
standoff),  and  364.0  m/s  at  5.09  ms  (2.00-m  standoff).  For  comparison,  the  curves  for  the 
thin  rectangular  and  trapezoidal  barricades  for  a  3.05-m  sta,ndoff  are  included  in  Figme  33 
with  the  same  legend  labels  as  were  used  for  their  momentum  plots  in  Figure  32.  For  the 
3  05-m  standoff  the  thin  rectangular  barricade  has  a  peak  velocity  of  476.6  m/s  and  the 
trapezoidal  barricade  has  a  peak  velocity  of  178.2  m/s.  The  X-direction  barricade  velocities 
become  less  meaningful  at  late  time  as  far  as  the  acceptor  stack  is  concerned  because  of 
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Table  2.  Barricade  Peak  X-Direction  Bulk  Motion  Parameters. 


Computation 

Number 

Massive 

Trape¬ 

zoidal 

980505 

Thin 
Rectan¬ 
gular 
980825  , 

Thick 

Rectan¬ 

gular 

980918 

Thick 

Rectan¬ 

gular 

980923 

Thick 

Rectan¬ 

gular 

980924 

Standoff  (m) 

3.048 

3.048  ^ 

3.048 

2.50 

2.00" 

Barricade 

Mass 

(kg/cm  of  depth) 

58.71 

28.61 

41.49 

41.49 

41.49 

Peak 

Momentum 
(Mg-m/s) 
Time  (ms) 

10.46 

10.07 

13.64 

6.560 

14.74 

8.334 

14.77 

6.559 

15.10 

5.089 

Peak 

Velocity 

(m/s) 

Time  (ins) 

178.2 

10.93 

476.6 

6.560 

355.3 

8.334 

355.9 

6.559 

364.0 

5.089 

Peak  Positive 
Acceleration 
(km/s/s) 

Time  (ms) 

143.4 

10.04 

631.5 

0.8141 

434.7 

0.8141 

453.3 

0.6947 

468.9 

0.5948 

Peak  Negative 
Acceleration 
(km/s/s) 

Time  (ms) 

-19.22 

32.80 

-808.7 

7.093 

-310.0 

9.494 

-317.8 

8.023 

-347.9 

6.883 

Peak 

Left-Surface 
Impulse 
(MN-s/m) 
Time  (ms) 

0.9016 

34.48 

1.185 

38.29 

1.241 

39.99 

1.308 

39.88 

1.396 
39.37  i 
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the  increasing  proportion  of  water  that  is 
acceptor  stack. 


in  the  air  above  the  plane  of  the  top  face  of  the 


Figure  33.  Water  Barricade  X-Direction  Velocity  Toward  the  Acceptor  Stack,  Computations 
980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and  98050j 
(Trapezoidal). 


The  momentum  curves  in  Figure  32,  and  hence  the  velocity  curves  in  Figure  33,  are  rel¬ 
atively  smooth  functions  with  respect  to  time.  The  velocities  were  piecewise  clifFerentiated 
with  respect  to  time,  using  the  difference  values  of  velocity  and  time  in  the  data  file,  to  pro¬ 
duce  the  curves  of  bulk  X-direction  acceleration  for  each  standoff  as  shown  in  Figure  34.  The 
peak  positive  accelerations  for  the  thick  rectangular  barricade  are  434.7  km/s^  at  0.81  ms 
(3.05-m  standoff),  453.3  km/s'^  at  0.69  ms  (2.50-m  standoff),  and  468.9  km/s^  at  0.59  ms 
(2.00-m  standoff).  For  comparison,  the  peak  positive  acceleration  of  the  thin  rectangular 
barricade  is  631.5  km/s^  at  0.814  ms.  For  the  trapezoidal  barricade  it  is  143.4  km/s  at 
10.0  ms.  These  plots  are  also  shown  in  Figure  34.  The  peak  negative  accelerations  (ne.,  de¬ 
celerations)  caused  by  the  barricade,  striking  the  acceptor  stack  are  particularly  informative, 
especially  when  viewed  in  conjunction  with  the  acceptor  stack  loading  that  is  presented  in  a 
following  section  of  this  report.  The  negative  accelerations  for  the  thick  rectangular  barricade 
are  minus  310  0  km/s‘'^  at  9.49  ms  (3.05-m  standoff),  minus  317.8  km/s^  at  8.02  ms  (2.50-m 
standoff),  and  minus  347.9  km/s^  at  6.88  ms  (2.00-m  sUndoff).  For  a^3.05-m  standoff,  the 
thin  rectangular  barricade  has  a  peak  negative  acceleration  of  minus  808.7  km/s  at  7.09  ms. 
The  trapezoidal  barricade  has  a  peak  negative  acceleration  of  minus  19.2  km/s  at  32.8  nas. 
Because  the  full  simulation  time  is  displayed  on  the  abscissa,  the  initial  accelerations  of  the 
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barricade  for  each  standoff  for  the  thick  rectangular  barricade  appear  to  nearly  overlay  one 
another.  Figure  35  shows  a  temporally  expanded  plot  of  the  first  12.0  ms  of  the  X-direction 
acceleration  of  the  barricade  for  each  computation.  The  initial  accelerations  for  the  thick 
rectangular  barricade  occur  in  a  direct  sequence  based  on  standoff  distance,  with  the  miti^ 
accelerations  for  the  thin  rectangular  and  the  trapezoidal  barricades  at^  a  3.05-m  stando 
beginning  at  about  the  same  time  as  that  for  the  thick  rectangular  barricade  at  that  same 
standoff.  The  deceleration  sequence  for  the  thick  rectangular  barricade  also  occurs  in  direct 
correspondence  to  the  standoff  distance.  No  meaningful  deceleration  of  the  trapezoidal  bar¬ 
ricade  occurs  during  the  first  12.0  ms,  but  the  thin  rectangular  barricade  clearly  shows  a 
great  deceleration  at  7.09  ms. 


Figure  34.  Water  Barricade  X-Direction  Acceleration  Toward  the  Acceptor  Stack,  Compu¬ 
tations  980918  Through  980924  (Thick  Rectangulai-),  Plus  980825  (Thin  Rectangular)  and 

980505  (Trapezoidal). 


Figure  36  shows  the  total  X-direction  impulse  per  meter  depth  on  the  left  surface  of 
the  barricade.  This  was  computed  by  integrating  the  overpressure  over  space  and  time 
using  30  tracer  particles  that  were  placed  along  the  left  surface  of  the  barricade  at  time 
zero  the  grid  generation  time;  The  overpressure  is  the  absolute  pressure  minus  the  ambient 
atmospheric  pressure.  The  tracer  particles  u^re  allowed  to  move  freely  with  the  flow  in 
the  grid  4s  the  simulated  time  in  the  computations  progressed,  the  left  surface  ol  the 
barricade  became  increasingly  distorted  to  the  point  that  it  was  no  longer  clearly  definable 
as  a  simple  surface.  Correspondingly,  the  impulse  integral  itself  probably  lost  meaning  after 
about  10  ms.  Essentially  all  of  the  impulse  from  the  detonation  of  the  donor  stack  is  delivered 
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p’igure  35.  Water  Barricade  Initial  X-Direction  Acceleration  Toward  the  Acceptor  Stack, 
Computations  980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangulai-) 
and  980505  (Trapezoidal). 


to  the  barricade  in  the  first  few  milliseconds.  There  is  a  moderate  inverse  functional  relation 
in  impulse  delivered  to  the  thick  rectangular  barricade  with  respect  to  standofi  distance. 
There  is  less  impulse  deli  vered  to  the  thin  rectangular  barricade  than  to  the  thick  rectangular 
barricade  at  a  3.05-m  standoff  distance.  Even  less  impulse  is  delivered  to  the  trapezoidal 
barricade  than  to  the  thick  rectangular  barricade  at  a  3.0.5-m  standoff  distance.  The  peak 
values  for  the  thick  rectangular  barricade  are  1.241  MN-s/m  at  40.0  ms  (3.0.5-m  standoff), 
1.308  MN-s/m  at  39.9  ms  (2.50-m  standoff),  and  1.396  MN-s/:rn  at  39.4  ms  (2.00-m  standoff). 
Because  of  the  surface  distortion  just  discussed,  the  times  of  these  peaks  are  not  particularly 
important  and  are  included  only  for  completeness.  This  equates  to  a  direct  ratio  of  peak 
impulse  for  the  thick  rectangular  barricade  of  1.125  for  an  inverse  ratio  in  relative  stancloff 
distance  of  1.524  for  the  standoff  range  of  3.05  m  to  2.00  m.  For  comparison,  the  peak 
impulse  for  the  thin  rectangular  barricade  is  1.185  MN-s/m  at  38.3  ms  (3.05-m  standoff), 
and  0.9016  MN-s/m  at  34.5  ms  for  the  trapezoidal  barricade  at  that  same  standoff. 

The  velocity  data  are  used  to  compute  the  bulk  translation  of  the  barricade  vemus 
time,  which  is  shown  in  Figure  37.  The  curves  for  the  thick  rectangular  barricade  are  most 
meaningful  through  about  6  to  7  ms.  At  about  this  time  the  barricade  impacts  the  acceptor 
stack.  The  peak  distances  for  the  thin  and  thick  rectangular  barricades  are  reduced  because 
of  the  net  rebound  of  parts  of  the  barricades  from  the  acceptor  stack  left  surface.  The 
curve  for  the  displacement  of  the  trapezoidal  barricade  in  Computation  980505  shows  less 
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Fieure  36  Water  Barricade  Left  Surface  Total  X-Direction  Impulse  per  Meter  Depth,  Com¬ 
putations  980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and 

980505  (Trapezoidal). 


displacement  through  about  30  ms  and  the  greatest  displacement  at  40  ms  because  there  is 
very  little  rebound. 

Figure  38  shows  the  functional  relations  of  the  peak  (at  diflerent  times)  values  of  the 
several  parameters  just  described  for  the  thick  rectangular  barricade  in  the  preceding  figures 
The  abscissa  shows  the  dimensional  standoff  distance.  The  ordinate  shows  the  normalizec 
direct  ratio  of  parameters,  the  value  of  a  given  parameter  at  a  given  standoff  divided  by  the 
corresponding  value  for  the  3.048-m  standoff.  Figure  39  shows  the  same  ordinate  data  as  m 
Figure  38  but  with  the  abscissa,  showing  the  normalized  inverse  standoff  ratio,  computed 
as  3  048  m  divided  by  each  successive  standoff.  Thus,  the  value  for  the  3.048-m  standoff 
itself  is  1  0,  and  the  value  for  the  2.00-m  standoff  is  1.524.  The  ordinate  is  the  same  as  for 
Figure  38  except  for  its  scaling.  The  abscissa,  and  ordinate  scales  are  forced  to  be  equal  so 
that  any  AX/ AT  =  ±1.0  relationship  would  show  as  a  ±4.5-degree  straight  line.  The  figures 
show  a  wea:k  fimctional  relation  of  all  of  these  normalized  parameters  with  both  dimensional 
and  normalized  standoff.  This  shows  that,  as  far  as  these  parameters  for  this  simplified 
barricade  are  concerned,  there  is  only  a.  minor  penalty  in  barricade  whole-bodj'  dynamics 
incurred  by  moving  the  barricade  closer  to  the  donor  stack  to  a  nominal  2-m  from  anommal 
.3-m  standoff.  The  apparent  lack  of  a  visible  black  curve  for  the  peak  momentum  m  both 
figures  is  caused  by  the  precise  overlay  by  the  red  curve  for  peak  velocity  because  of  their 

direct  scaling  by  mass. 
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Figure  37  Water  Barricade  X-Direction  Distance  Moved  Toward  tlie  Acceptor  Stack,  Com¬ 
putations  980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and 
980505  (Trapezoidal). 


Figure  38.  Normalized  (Direct  Ratio)  Thick  Rectangular  Barricade  Parameters  Versus 
Standoff  Distance,  Computations  980918  Through  980924. 
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Figure  39:  Normalized  (Direct  Ratio)  Thick  Rectangular  Barricade  Parameters  Versus  Nor¬ 
malized  (Indirect  Ratio)  Standoff  Distance,  Computations  980918  Through  980924. 


3.3.  Acceptor  Stack  Dynamics 

Figure  40  shows  the  bulk  momentum  per  centimeter  depth  of  the  acceptor  stack  in  the  X 
direction  for  the  three  computations  for  the  thick  rectangular  barricade,  plus  Cornputations 
980825  for  the  thin  rectangular  barricade^  and  980505  for  the  massive  trapezoidal  barNcade, 
both  of  which  were  for  a  3.05-m  standoff.  Positive  momentum  is  defined  in  the  positive  X 
direction  as  before.  There  is  only  a  minimal  increase  in  the  momentum  of  the  acceptor  stack 
caused  by  the  air  shock  for  the  three  thick  rectangular  barricade  computations.  .After  about 
5  ms,  all  three  computations  show  a  very  rapid  increase  in  momentum  caused  by  the  impact 
of  the  water  bairicade.  The  acceptor  stack  momentum  for  the  thin  rectangular  barricade 
from  the  3.05-m  standoff  computation  begins  a  rapid  increase  about  2  ms  sooner  than  the 
corresponding  curve  for  the  thick  rectangular  barricade  at  the  same  standoff  and  reaches 
the  greatest  value  of  11.90  Mg-m/s  for  all  curves  by  40  ms.  The  acceptor  stack  momentum 
from  the  3.05-m  standoff  computation  for  the  massive  trapezoidal  water  barricade  shows  the 
three-stage  secpience  increase  in  momentum  described  previously^  that  is  caused  by  the  air 
shock,  the  impact  of  the  water  wave  on  the  top-left  surface  of  the  acceptor  stack,  and  then 
the  impact  of  the  lower  sections  of  the  barricade  on  the  acceptor  stack.  The  momentum  of 
the  acceptor  stack  at  a  nominal  time  of  40  ms  for  Computation  980918  (thick  rectangular 
barricade.  3.05-m  standoff,  28.61  kg/cra  of  depth)  is  11.11  Mg-m/s,  which  is  2.80  times 
o-reater  than  the  corresponding  value  of  3.962  Mg-m/s  for  Gomputation  980505  (trapezoidal 
barricade,  3.05-m  standoff,  5S.T1  kg/cm  of  depth).  Theffnal  values  of  the  momentum  of  the 


acceptor  stack  for  the  thick  rectangular  barricade  computations  at  a  nominal  time  of  40  ms 
are  11.35  Mg-m/s  for  the  2.50-m  standoff  and  11.73  Mg-m/s  for  the  2.00-m  standoff.  Table  3 
contains  a  summary  of  several  X-direction  parameters  that  describe  some  of  the  bulk  motion 
of  the  acceptor  stack  for  the  various  computations.  First  among  those  parameters,  after 
the  computation  numbers  and  standoff  distances,  are  the  peak  X-direction  bulk  momentum 
values  for  the  acceptor  stack,  along  with  their  respective  times  of  occurrence,  listed  with  more 
significant  figures  than  were  typically  used  in  the  text  for  completeness.  In  order  to  facilitate 
comparisons,  the  first  column  of  numbers  is  for  the  acceptor  stack  in  Computation  980505 
for  the  massive  trapezoidal  barricade  at  a  3.05-m  standoff,  followed  by  data  for  Computation 
■  980825  for  the  thin  rectangular  barricade  at  a  3.05-m  standoff,  and  then  followed  by  columns 
of  data  for  the  three  thick  rectangular  barricade  computations.  The  rest  of  the  parameters 
in  the  table  are  discussed  in  the  following  paragraphs. 


Figure  40.  Acceptor  Stack  X-Direction  Momentum,  Computations  980918  Through  980924 
(Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and  980505  (Trapezoidal). 

The  corrected  mass  of  the  acceptor  stack  was  used  to  compute  the  bulk  X-direction 
velocity  from  the  momentum  of  the  acceptor  stack.  The  results  are  shown  in  Figure  41.  The 
curves  show'  the  same  timing  and  differentiation  as  those  for  the  acceptor  stack  momentum. 
The  velocities  at  40  ms  for  the  thick  rectangular  barricade  are  93.66  m/s  (3.05-m  standoff), 
95.68  m/s  (2.50-m  standoff),  and  98.93  m/s  (2.00  m  standoff).  The  velocities  at  40  ms  at  a 
3.0.5-m  standoff  for  the  thin  rectangular  water  barricade  and  the  massive  trapezoidal  w'ater 
barricade,  also  shown  in  Figure  41,  are  100.3  m/s  and  33.4  m/s,  respectively. 


Table  3.  Acceptor  Stack  Peak  X-Directioii  Bulk  Motion  Parameters. 


Computation 

Number 

Massive 

Trape¬ 

zoidal 

980505 

Thin 
Rectan¬ 
gular 
980825  ' 

Thick 

Rectan¬ 

gular 

980918 

Thick 

Rectan¬ 

gular 

980923 

Thick 

Rectan¬ 

gular 

980924 

Standoff  (m) 

3.048 

3.048 

3.048 

2.50 

2.00 

Peak 

Momentum 
(Mg-m/s) 
Time  (ras) 

3.962 

40.00 

11.90 
39.87  ■ 

11.11 

39.62 

11.35 

39.71 

11.73  ' 

39.69 

Peak 
Velocity 
(m/s) 
Time  (rns) 

33.40  1 
40.00  ' 

100.3 

39.87 

93.66 

39.62 

95.68 

39.71 

98.93 

39.69 

Peak 

Acceleration 
(km/s/s) 
Time  (ms) 

9.277 
32.80  : 

197.5 

7.093 

109.8  1 
9.494 

112.8 

8.023 

123.6 

6.883 

Peak 

Left-Surface 
Impulse 
(MN-s/m) 
Time  (ms) 

0.3725 

39.99 

1.138 

39.99 

1.063 

39.99 

1.098 

39.99 

1.144 

39.99 

Distance 

Traveled 

(m) 

Time  (ms) 

0.4065 

39.99 

3.237 

39.99 

2.805 

39.99 

3.014 

39.99 

3.246  : 
39.99 
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Figure  41.  Acceptor  Stack  X-Direction  Velocity,  Computations  980918  Through  980924 
(Thick  RectangulaT),  Plus  980825  (Thin  Rectangular)  and  980505  (Trapezoidal). 
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As  was  done  for  the  barricade,  the  acceptor  stack  velocity  for  each  standoff  was  piecewise 
differentiated  with  respect  to  time  to  compute  the  bulk  X-direction  acceleration  of  the  ac¬ 
ceptor  stack.  The  acceleration  curves  for  the  three  thick  rectangular  barricade  computations 
and  Computations  980505  and  980825  are  shown  in  Figure  42.  Each  individual  curve  for 
the  thick  rectangular  barricade  shows  a  large  spike  in  acceleration  in  the  order  of  increasing 
standoff  distance.  All  occur  before  10  ms.  The  curves  for  the  2.50-m  and  2. 00-m  standoffs 
show  some  leading  structure  in  the  acceleration  curves  prior  to  the  main  acceleration  spike, 
most  likely  because  of  the  leading  impact  of  the  bottom  section  of  the  barricade  against  the 
bottom  of  the  left  face  of  the  acceptor  stack,  followed  by  a  progressive  impact  of  the  rest 
of  the  barricade  moving  up  the  acceptor  stack  left  face.  The  curve  for  the  thin  rectangular 
barricade  at  a  3.05-m  standoff  shows  a  single,  large  spike  in  acceleration  at  an  early  time,  as 
does  the  curve  for  the  thick  rectangular  barricade  at  that  same  standoff.  The^peak  accelera¬ 
tions  for  the  two  different  rectangular  barricades  are  far  greater  than  the  9.28  km/s‘^  for  the 
trapezoidal  barricade  at  a  3.05-m  standoff  that  occurs  at  32.8  ms. 


Figure  42.  Acceptor  Stack  X-Direction  Acceleration,  Computations  980918  Through  980924 
(Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and  980505  (Trapezoidal). 

The  acceptor  stack  was  modeled  as  a  solid  iron  rectangle  so  that  the  most  reliable 
loading  possible  could  be  computed  for  its  left  face.  Thirty  tracer  particles  were  uniformly 
spaced  along  the  left  face,  top  to  bottom',  of  the  acceptor  stack.  They  were  constrained  from 
mo\-ing  in  either  the  X  or  Y  direction  so  that  the  ensuing  hydrodynamic  flows  would  not 
sweep  them  off  the  face  of  the  acceptor  stack  or  reposition  them  horizontally  or  vertically. 
The  overpressure  histories  were  integrated  over  space  and  time  to  compute  the  total  X- 
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direction  impulse  per  meter  depth  versus  time  for  each  standoff.  These  curves  are  shown  in 
Figure  43.  The  acceptor  stack  shows  that  the  thick  rectangular  water  barricade  provided  a 
very  efficient  delivery  of  its  left-face  impulse  from  the  blast  loading  caused  by  the  detonation 
of  the  donor  stack  (see  Figure  36)  to  the  acceptor  stack  at  its  right  face.  The  acceptor  stack 
left-face  values  are  nearly  equal  to  the  respective  barricade  left-face  values.  The  final  impulse 
values  per  meter  depth  on  the  acceptor  stack  left  face  are  1.063  MN-s/m  (3.05-m  standoff), 
1.098  MN-s/m  (2.50-m  standoff),  and  1.144  MN-s/m  (2.00-m  standoff).  For  comparison,  the 
final  impulse  values  on  the  left  face  of  the  thick  rectangular  water  barricade  (see  Table  2) 
are  1.241  MN-s/m  (3.05-m  standoff),  1.308  MN-s/m  (2.50-m  standoff),  and  1.396  MN-s/m 
(2.00-m  standoff).  The  term  “impulse-transfer  efficiency”  is  herein  defined  as  the  impulse 
delivered  to  the  left  face  of  the  acceptor  stack  divided  by  the  impulse  delivered  to  the  left  face 
of  the  barricade.  The  impulse-transfer  characteristics  of  all  of  the  water  barricades  studied 
to  date  are  summarized  in  Table  4.  Although  impulse  and  impulse-transfer  efficiency  as 
defined  and  used  herein  are  not  the  only  valid  indicators  of  the  efficacy  of  a  given  barricade 
design,  they  are  important.  The  data  listed  in  Table  4  are  grouped  by  barricade  design  in 
ascending  order  of  impulse-transfer  efficiency,  moving  from  top  to  bottom  in  the  table.  A 
low  value  of  impulse  on  the  acceptor  stack  is  desirable,  as  is  a  low  value  of  impulse-transfer 
efficiency.  The  thin  rectangular  barricade  showed  the  worst  performance  by  this  measure, 
having  the  highest  values  of  impulse-transfer  efficiency.  The  rate  of  delivery  of  impulse  is  also 
very  important  because  of  its  direct  relation  to  peak  loads  and  accelerations.  Two  simple 
but  relatively  effective  methods  for  showing  the  peak  impulse  on  the  left  face  of  the  acceptoi 
stack  versus  a  geometrically  scaled  barricade  mass  were  found,  which  help  to  tie  the  results 
of  all  of  the  water  barricade  computations  together.  The  first  method,  Scaling  Method  1,  is 
shown  in  Figure  44.  The  water  barricade  mass,  in  units  of  kg/cm  of  depth,  was  multiplied 
by  the  square  of  the  secant  of  the  angle,  0,  to  the  vertical  of  the  front  and  back  faces,  once 
for  the  slope  of  the  front  face  and  again  for  the  slope  of  the  back  face  for  a  compounded 
factor  equal  to  the  fourth  power  of  the  secant  of  0.  This  produced  a  scaled  barricade  mass, 
Ml,  from  the  mass  in  kg  per  cm  of  depth,  M,  and  the  face  angle  in  degrees,  0,  shown  in 


equation  1. 

Ml  =  M  (sec  (0))^ 


(1) 


For  the  trapezoidal  water  barricade,  0  is  equal  to  30  degrees  for  both  faces.  This  resulted 
in  the  three-tiered  set  of  lines  shown  in  Figure  44,  one  line  for  each  standoff  distance.  A 
second  method.  Scaling  Method  2,  is  shown  in  Figure  45,  in  which  the  scaled  bamcade  mass 
in  Figure  44  is  further  multiplied  by  the  cube  root  of  the  standoff  distance,  S',  in  meters  to 
form  a  new  scaled  barricade  mass,  M2.  This  is  shown  in  equation  2. 


(2) 


The  data  are  shown  with  symbols  to  make  it  easier  to  discuss  the  two  differentiated  sets 
that  are  indicated  in  this  figure.  The  first  set,  represented  by  those  symbols  for  which  the 
scaled  abscissa  value  is  less  than  75,  is  comprised  of  points  only  for  the  rectangular  water 
barricade.  The  second  set,  represented  by  those  symbols  for  which  the  scaled  abscissa^ value 
is  greater  than  125,  is  comprised  of  points  only  for  the  trapezoidal  water  barricade.  These 
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two  sets  fall  almost  on  the  same  line.  These  relatively  simple  and  reasonably  successful 
attempts  at  scaling  the  impulse  data  were  not  pursued  any  further  because  they  are  based 
only  on  a  limited  set  of  data  for  water.  No  further  scaling  is  proposed  at  this  time  until 
additional  computations  can  be  made  for  sand  barricades,  at  which  time  factors  for  density 
and  perhaps  other  shapes  should  be  considered. 


Figure  4.3.  Acceptor  Stack  X-Direction  Total  Impulse  per  Meter  Depth,  Computations 
980918  Through  980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and  980o0o 

(Trapezoidal). 


Figure  46  shows  the  distance  that  the  acceptor  stack  moves  as  a  result  of  the  blast  and 
impact  loading  bv  40.0  ms.  The  range  for  Computations  980918  through  980924  for  the 
thick  rectangular-  barricade  is  from  2.80  m  (3.05-m  standoff)  to  3.25  m  (2.00-m  standoff). 
Computation  980825  for  the  thin  rectangular  barricade  shows  a  movement  of  the  acceptor 
stack  of  3.24  m  at  the  3.05-m  standoff.  In  Computation  980505  for  the  trapezoidal  barricade, 
the  acceptor  stack  moves  0.41  m  in  40.0  ms,  12.6  percent  of  that  for  the  acceptor  stack  in 
Computation  980825,  and  14.5  percent  of  that  for  Computation  980918. 

Figure  47  shows  the  functional  relations  of  the  peak  (at  different  times)  and  final  (at 
40  0  ms)  values  of  the  several  parameters  versus  standoff  distance  that  were  just  described 
for  the  acceptor  stack  with  the  thick  rectangular  barricade  in  the  preceding  figures.  The 
ordinate  parameters  are  normalized  in  the  same  way  as  was  done  m  Figure  38:  the  c  irect 
ratio  of  the  respective  parameters  relative  to  the  values  for  the  3.05-m  standoff.  1  he  black 
momentum  curve  does  not  show  because  it  is  overlaid  by  the  red  velocity  curve.  This  is 
because  of  the  simple  scaling  of  the  two  curves  by  mass.  Figure  48  shows  the  same  ordinate 
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Table  4.  Impulse-Transfer  Efficiency  of  All  Water  Barricades. 


Barricade 

and 

Standoff 

Peak  Barricade 
Left-Surface 
Impulse 
(MN-s/m) 

Peak  Acceptor  Stack 
Left  Surface 
Impulse 
(MN-s/m) 

Impulse- 

Transfer 

Efficiency 

(-) 

Massive  Trapezoidal 
at  3.048-m  Standoff 

0.9016 

0.3725 

0.4132 

Massive  Trapezoidal 
at  2.75-m  Standoff 

0.9337 

0.4044 

0.4331 

Massive  Trapezoidal 
at  2.50-m  Standoff 

0.9381 

0.4195 

0.4472 

Massive  Trapezoidal 
at  2.25-m  Standoff 

0.9502 

0.4279 

0.4503 

Massive  Trapezoidal 
at  2.00-m  Standoff 

0.9630 

0.4945 

0..5135 

Thick  Rectangular 
at  3.048-m  Standoff 

1.241 

1.063 

0.8.566 

Thick  Rectangular  n 

at  2.50-m  Standoff 

1.308 

1.098 

0.8395 

Thick  Rectangular 
at  2.00-m  Standoff 

1.396 

1.144 

Thin  Rectangular 
at  3.048-m  Standoff 

1.185 

1.1.38 

0.9603 

Thin  Rectangular 
at  2.50-m  Standoff 

1.283 

1.154 

0.8995 

Thin  Rectangular 
at  2.00-m  Standoff 

1.360 

1.248 

0.9177 
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Scaled  Barricade  Mass.  M,  {kg/cm)(-) 

Figure  44.  Acceptor  Stack  Peak  X-Direction  Total  Impulse  per  Meter  Depth  Versus  Barn 
cade  Mass,  All  Water  Barricade  Computations,  Scaling  Method  1. 


Scaled  Barricade  Mass,  iMg  {kg/cm){-)(m) 


Figure  45.  Acceptor  Stack  Peak  X-Direction  Total  Impulse  per  Meter  Depth  Wrsus  Ba.rri 
cade  Mass,  All  Water  Barricade  Computations,  Scaling  Method  2. 
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Figure  46.  Acceptor  Stack  X-Direction  Distance  Moved,  Computations  980918  Through 
980924  (Thick  Rectangular),  Plus  980825  (Thin  Rectangular)  and  980505  (Trapezoidal). 

data,  plotted  against  the  inverse  normalized  standoff,  with  that  normalization  done  in  the 
same  way  as  for  P’igure  S9,  including  a  forced  scaling  of  the  ordinate  and  abscissa  to  have  each 
cover  the  same  range  of  normalized  data  over  the  same  a,xis  length  to  facilitate  comparison 
with  Figure  39.  Like  the  barricade,  the  normalized  acceptor  stack  ordinate  parameters  are 
all  relatively  w'eak  direct  functions  of  inverse  standoff  ratio. 

The  simple  standoff  distance,  measured  from  the  stack  base  to  the  barricade  base,  may 
not  be  the  only  meaningful  distance  to  consider  when  examining  these  parameters  for  the 
acceptor  stack.  The  total  distance  between  the  right  face  of  the  donor  stack  and  the  left  face 
of  the  acceptor  stack,  equal  to  twice  the  standoff  plus  the  base  width  of  the  bamcade,  may  be 
an  informative  parameter  to  use.  For  convenience,  this  distance  is  hereinafter  referred  to  as 
‘'face  separation.”  Figure  49  shows  the  ordinate  values  from  Figures  47  and  48  plotted  against 
an  abscissa  showing  the  face  separation.  Figure  50  shows  the  same  ordinate  data  plotted 
against  the  inverse  normalized  face  separation.  The  normalizing  value  in  the  numerator  was 
the  face  separation  for  the  3.05  m  standoff.  One-to-one  scaling  for  the  normalized  abscissa 
and  ordinate  was  forced  in  this  .figure.  All  parameters  show  a  weak  correlation  with  the 
inverse  normalized  face  separation. 
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Figure  47.  Normalized  (Direct  Ratio)  Acceptor  Stack  Parameters  Versus  Standoff  Distance, 
Computations  9S091S  Through  980924  (Thick  Rectangular  Barricade). 


Figure  48.  Normalized  (Direct  Ratio)  Acceptor  Stack  Parameters  Versus  Normalized  (Inverse 
Ratio)  Standoff  Distance,  Computations  980918  Through  980924. 
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Figure  49.  Normalized  {Direct  Ratio)  Acceptor  Stack  Parameters  Versus  Face  Separation, 
Computations  980918  Througli  980924. 


Figure  50.  Normalized  (Direct  Ratio)  Acceptor  Stack  Parameters  Versus  Normalized  (Inverse 
Ratio)  Face  Separation,  Computations  980918  Through  980924. 
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3.4.  Acceptor  Stack  Left  Surface  Pressures 


The  pressures  on  the  surface  of  the  acceptor  stack  during  this  type  of  event  are  of  great 
interest  Data  from  the  30  tracers  that  were  placed  uniformly  along  the  left  surface  of 
the  acceptor  stack  were  processed  to  present  a  comprehensive  summary  of  the  overpressure 
history  on  that  surface  for  each  computation.  An  area-weighted  averap  overpressure  was 
computed  using  all  of  the  30  individual  tracer  pressures  at  each  point  in  time.  The  ma.xiimra 
overpressure  for  any  of  the  tracers  at  a  given  time  was  identified,  as  was  the  minimum.  The 
results  from  the  three  computations  for  the  thick  rectangular  barricade  are  presented  here  for 
each  standoff  distance.  Figure  51  shows  the  average,  maximum,  and  minimum  overpressures 
versus  time  on  the  left  surface  of  the  acceptor  stack  for  a  standoff  of  3.05  m  for  Computation 
980918  It  shows  a  peak  average  overpressure  of  0.524  GPa  at  9.49  ms  and  a  peak  overpressure 
for  an  individual  point  of  2.49  GPa,  also  at  9.49  ms.  Similarly,  Figure  52  shows  the  average, 
maximum,  and  minimum  overpressures  versus  time  on  the  left  surface  of  the  acceptor  stack 
for  a  standoff  of  2.50  m  for  Computation  980923.  It  shows  a  peak  average  overpressure 
of  0  536  GPa  at  8.02  ms  and  a  peak  overpressure  for  an  individual  point  of  2.00  GPa  at 
7.75  ms.  Figure  53  shows  the  average,  maximum,  and  minimum  overpressures  vemus  time 
on  the  left  surface  of  the  acceptor  stack  for  a  standoff  of  2.00  m  for  Computation  980924.  It 
shows  a  peak  average  overpressure  of  0.580  GPa  at  6.84  ms  and  a  peaJ<  overpressure  for  an 
individual  point  of  1.76  GPa  at  6.71  ms.  There  is  an  interesting,  somewhat  contradictory 
trend  of  increasing  average  overpressure  on  the  acceptor  stack  left  surface  with  decreasing 
standoff  distance  coincident  with  a  trend  of  decreasing  peak  overpressure  for  an  individual 
point  with  decreasing  standoff  distance.  It  appears  to  be  caused  by  a  trade-off  in  the  systein 
dynamics  in  the  bulk  acceleration  of  the  barricade,  and  in  the  time  and  distance  required  to 
accelerate  it  before  it  strikes  the  acceptoi  stack. 


The  figures  that  follow  each  show  the  average  overpressure  and  the  peak  overpressure 
for  an  individual  point  for  all  water  barricade  computations  reported  to  date  for  the  massive 
trapezoidal  water  barricadel  ^^d  the  thin  rectangular  water  barricade^  at  matching  standoff 
distances  to  the  values  from  Computations  980918  through  980924  just  shown  in  Figures  ol, 
and  53  Figure  54  show's  the  average  and  individual  peak  overpressures  versus  time  on 
the  left  surface  of  the  acceptor  stack  for  a  standoff  of  3.05  m  for  Computations  980918, 
980825,  and  980505  that  are  for  the  thick  rectangular,  thin  rectangular,  and  trapezoida 
water  barricades,  respectively.  The  order,  from  highest  to  lowest,  of  the  peak  individual 
overpressure  is  thick  rectangular,  thin  rectangular,  and  then  trapezoidal.  The  order  for  the 
average  overpressure  is  thin  rectangular  (greatest),  thick  rectangular,  and  trapezoidal,  ig- 
ure  55  shows  the  average  and  individual  peak  overpressures  versus  time  on  the  left  surtace 
of  the  acceptor  stack  for  a  standoff  of  2.50  m  for  Computations  980923,  980826,  and  980oLl 
that  are  for  the  thick  rectangular,  thin  rectangular,  and  trapezoidal  w'ater  barricades  re¬ 
spectively.  The  order,  from  highest  to  lowest,  of  the  peak  individual  overpressure  is  thick 
rectanguiar  (bv  a  small  margin),  thin  rectangular,  and  then  trapezoidal.  The  order  for  the 
average  overpressure  is  thin  rectangular  (greatest),  thick  rectangular,  and  trapezoidal.  Fig¬ 
ure  56  shows  the  average  and  individual  peak  overpressures  versus  time  on  the  left  surface 
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Figure  51.  Acceptor  Stack  Left  Surface  Overpressure,  3.05-m  Standoff,  Computation  980918. 


Figure  52.  Acceptor  Stack  Left  Surface  Overpressure,  2.50-m  Standoff,  Computation  980923. 


Figure  53.  Acceptor  Stack  Left  Surface  Overpressure,  2.00-m  Standoff,  Computation  980924. 

of  the  acceptor  stack  for  a  sta.ndoff  of  2.50  m  for  Computations  980924,  980827,  and  980610 
that  are  for  tlie  thick  rectangular,  thin  rectangular,  and  trapezoidal  water  barncades,  re¬ 
spectively  The  order,  from  highest  to  lowest,  of  the  peak  individual  overpressure  is  thin 
rectanguiar,  thick  rectangular,  and  then  trapezoidal.  The  order  for  the  average  overpressure 
is  thin  rectangular  (greatest),  thick  rectangular,  and  trapezoidal.  For  all  plots,  the  mas¬ 
sive  trapezoidal  barricade  always  showed  the  smallest  values  by  far  of  both  peak  individual 
overpressure  and  average  overpressure  on  the  acceptor  stack  left  face.  All  plots  showed  an 
essentially  negligible  loading  from  the  air  shock  at  earfy  time.  All  computations  for  bo  h 
the  thin  and  the  thick  rectangular  wa.ter  barricades  show  peak  vaJues  of  overpressure  on  the 
left  surface  of  the  acceptor  stack  in  the  range  of  1.75  CiPa  (17.5  kbar)  or  higher  These  are 
high  enough  pressures  to  represent  a  threat  of  inducing  a  chemical  rearfion  in  the  acceptor 
sta.ck  if  they  are  efficiently  transmitted  through  packaging  and/or  casings  to  the  eiiergetic 
loads  of  the  munitions.  The  report  by  Liddiard  and  Forbesl-^  stated,  for  example,  that  the 
underwater  sensitivity  test  (UST)  showed  that  "...compression  by  a  3  or  4  kbar  shock  is 
of  itself  a  sufficient  external  stimulus  to  start  chemical  reaction  in  a  heterogeneous  solid 
explosive  such  as  pentolite...”  and  "...UST  burning  occurs  at  peak  stresses  of  4  to  12  kbar 
in  the  explosives...”  A  cautionary  note:  a  simple  shock  stimulus  is  not  the  only  initiating 

mechanism  for  an  explosive. 
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Figure  54.  Acceptor  Staclc  Left  Surface  Overpressure,  3.05-m  Staucloff,  Computations  980918 
(Thick  Rectangular),  980825  (Thin  Rectangular)  and  980505  (Trapezoidal). 


Figure  55.  Acceptor  Stack  Left  Surface  Oveiijressure,  2.50-in  Standoff,  Computations  980923 
(Thick  Rectangular),  980826  (Thin  Rectangular)  and  980521  (Trapezoidal).. 


67 


Figure  56.  Acceptor  Stack  Left  Surface  Overpressure,  2.00-m  Standoff,  Computations  980924 
(Thick  Rectangular),  980827  (Thin  Rectangular)  and  980610  (Trapezoidal). 


4.  CONCLUSION 

The  coupled  computations  discussed  herein  modeled  a  simplified,  uncased,  rectangulai 
explosive  chai-ge  representing  a  nominal  munitions  stack  containing  4,000  kg  of  Comp-B 
undergoing  a  complete,  high-order  detonation  with  the  initiation  point  at  its  center  iNo  mu¬ 
nitions  casings  or  packing  materials  (and  their  resulting  fragments)  were  included.  The  only 
harricade  design  that  was  used  was  a  solid,  water-only  1.70-m-thick  rectangle.  These  compu¬ 
tations  were  compared  with  computations  for  a  thin,  1.17-m  rectangular  water  barricade  and 
a  massive  trapezoidal  water  barricade  at  the  same  standoff  distances.  Those  computations 
were  reported  previously.  A  geometrically  simplified  2-D  Cartesian  coordinates  system  with 
the  same  finite-difference  grid  was  used  throughout  the  computations.  This  eliminated  three- 
dimensional  divergence  effects  that  could  reduce  loadings  considerably.  The  only  parameter 
that  was  varied  was  the  standoff  distance. 

The  computations  for  the  thick  rectangular  water  barricade  demonstrate  a  relatively 
wea.k  inverse  functional  relationship  between  normalized  values  of  the  standoff  distance  and 
the  loading  on  and  whole-body  response  of  the  barricade.  Similar  results  for  both  stanc  o 
and  face  separation  were  found  for  the  loading  on  and  whole-body  response  for  the  acceptor 
stack.  The  impact  loading  on  the  acceptor  stack  by  the  thick  rectangular  water  barricade  is 
much  more  severe  than  that  reported  earlier^  for  the  massive  trapezoidal  water  barricade 
but  somewhat  less  severe  than  for  the  thin  rectangular  water  barricade.  Peak  individual 
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overpressures  on  the  acceptor  stack  for  the  two  different  rectangular  water  barricades  are 
high  enough,  1.75  GPa  or  greater,  to  represent  a  threat  of  initiating  a  chemical  reaction  in 
munitions  within  the  acceptor  stack.  Also,  neither  rectangular  barricade  was  effective  in 
keeping  explosive  products  from  the  donor  stack  away  from  the  acceptor  stack.  The  trape¬ 
zoidal  barricade  was  effective  in  that  regard.  A  relatively  simple  scaling  of  barricade  mass 
with  geometric  parameters  resulted  in  developing  a  single  functional  relationship  with  peak 
impulse  for  all  water  barricade  shapes  and  standoff  distances  considered  in  the  complete 
set  of  computations.  Because  the  donor  stack  was  represented  by  a  simple,  bare  explosive 
charge,  the  synergistic  effects  of  the  impact  of  large  numbers  of  high-speed  fragments  along 
with  the  barricade  impact  loading  were  not  addressed.  Additional  computational  studies  of 
sand-filled  barricades  were  in  progress  as  of  the  expiration  of  this  customer  project.  The 
possibility  of  extending  the  project  to  perform  the  sand  barricade  computations  and  com¬ 
putations  simulating  the  impact  of  water  or  sand  on  individual  munitions  is  being  explored. 
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avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to  organization, 
technical  content,  format,  etc.) _ _ 


Organization 

CURRENT  Name  E-mail  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  maU.) 
(DO  NOT  STAPLE) 


